Non-Abelian gauge potentials in graphene bilayers 
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We study the effect of spatial modulations in the interlayer hopping of graphene bilayers, such as 
those that arise upon shearing or twisting. We show that their single-particle physics, characterized 
by charge accumulation and recurrent formation of zero-energy bands as the pattern period L 
increases, is governed by a non-AbeUan gauge potential arising in the low-energy electronic theory 
due to the coupling between layers. We show that such gauge-type couphngs give rise to a potential 
that, for certain discrete values of L, spatially confines states at zero energy in particular regions 
of the Moiré patterns. We also draw the connection between the recurrence of the fiat zero-energy 
bands and the non-Abelian character of the potential. 



Introduction. — The discovery of graphene, the mate- 
rial made of a one-atom-thick carbon layer, has provided 
the realization of a system where the electrons have con- 
ical valence and conduction bands, therefore behaving as 
massless Dirac fermions [THS]- A remarkable feature of 
graphene is that deformations of its honeycomb lattice 
may produce a similar effect to that of gauge potentials 
in the low-energy Dirac theory [4] . Recently, it has been 
shown that the locai in-plane deformations induced by 
strain can be mimicked by an effective vector potential, 
which may give rise to the analogue of Landau levels in 
the deformed graphene sheet [S]. 

In this paper we show that the effect of modulations in 
the interlayer hopping of graphene bilayers can be rep- 
resented in general by a non-Abelian background gauge 
potential in the low-energy electronic theory, and that it 
is responsible for the zero energy charge density waves 
and the dispersionless minibands, predicted by theory, 
and recently measured [6]. The vector components of 
the potential take values in the space of SU(2) matrices, 
which correspond to rotations in the Hilbert space of the 
two layers. This kind of non-Abelian gauge fields [7j is 
relatively rare in a condensed-matter context [SHTO] , but 
it is quite relevant in subatomic physics, being responsi- 
ble for the interaction between matter fields. The proton 
and the neutron, for instance, compose an isospin SU(2) 
doublet. It was proposed long ago that an ideal experi- 
ment of scattering of these particles onto a non-Abelian 
flux line should lead to the transfer of protons into neu- 
trons and vice versa [TT]. In general, matter fields pick 
up a matrix-valued 'phase' in their propagation in a non- 
Abelian gauge fleld. Interference of such matrix-valued 
phase along two indistinguishable paths (as opposed to 
the conventional U(l) phase) leads to an intriguing non- 
Abelian generalization of the Aharonov-Bohm effect. In 
our context, this would manifest as coherent layer polar- 
ization induced by the interference of two SU(2) phases 
acquired along the two paths. 

Experimental realizations of non-Abelian gauge poten- 
tials have been proposed before in the study of ultracold 
atoms [IH [13 ■ Investigations have addressed in partic- 



ular the infìuence of the non-Abelian gauge potentials 
in the development of the Landau levels produced by a 
conventional magnetic field [l4 [ I15 j . However, the ques- 
tion of whether pure non-Abelian gauge fields may lead 
to a phenomenology similar to the magnetic localization 
of Landau states remains open. Our investigation sheds 
light on this question, showing that it is possible to de- 
velop a zero-energy level of spatially confined states as 
a consequence of the non-Abelian gauge potential, pro- 
vided that the fermion fields return to the originai inter- 
nai state around a closed path. 

The effective non-Abelian gauge potentials that arise in 
the bilayers have actually a genuine applied interest, since 
they induce periodic spatial confinement of electronic 
states. Indeed, we will see that the one-dimensional (ID) 
modulation of the interlayer tunneling leads to the con- 
finement of electronic states into narrow ID channels. We 
will also extend our approach to the description of twisted 
bilayers [T1H21], where the non-Abelian gauge potential 
turns out to confine low energy electrons into a trian- 
gular array of quantum dots. The problem of confine- 
ment of electronic states has particular relevance, given 
that scalar potential barriers are not effective to constrain 
the propagation of the electrons in graphene |25| , which 
makes the non-Abelian gauge potentials proposed in this 
paper an interesting alternative to the confinement and 
manipulation of electronic states in graphene devices. 

Model. — The simplest realizations of a non-Abelian 
gauge potential are found by means of a modulated mis- 
match in the relative position of the two lattices of a 
bilayer, obtained either by applying strain or shear in 
one of the layers or by relative rotation between the two 
layers. In both instances, the resulting mismatch pro- 
duces characteristic Moiré patterns, see Fig. [T] which 
refiect the spatial alternation between A^'-type stack- 
ing (perfect alignment of the atoms in the two layers) 
and AB'-type, i?A'-type (Bernal) stacking, where A{') 
and i?(') correspond to the two sublattices of the lower 
(upper) lattice. 

At energies e < 1 eV, the Moiré electron system is 
described by Dirac fermions on each layer, coupled by 
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FIG. 1. Moiré patterns of (a), top, sheared bilayer (showing 
the alternation between AA' , AB' and BA' stackings, and (b) 
twisted bilayer, where the hexagonal supercell and the difFer- 
ent types of stacking have been marked. (a), bottom, shows 
the eflective potential Vc{t{x) arising from the non-Abelian 
gauge potential A, together with a typical zero-energy state 
confined between the AB' and BA' regions, and a finite en- 
ergy state concentrated around AA' . 



a position dependent interlayer hopping amplitude. The 
Hamiltonian takes the forni [121 US] 
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where n± = —idx+dy^(Ax + iAy). The spatially modu- 
lated interlayer coupling functions V arise from the Moiré 
pattern formation, and the intra-layer Abelian gauge field 
±A describes the strains in each layer. These strains 
lead to Constant gauge fields in our case, A = AK/2. 
Note that, as discussed later, this model also describes 
twisted bilayers, in which the interlayer Dirac cone shift 
AK arises due to the relative twist between layers, not 
strains. Since A is anyhow uniform, we can gauge it 
away by a transformation U — exp((j/2)T3AK • r), where 
T3 is a Pauli matrix which operates on the layer in- 
dex. This transforms consequently the interlayer cou- 
plings into Vij{r) = Vij{r)e-'^^ '' . 

Low-energy theory of sheared bilayers. — We consider 
first the instance in which shear u^y is applied along 
the AB bonds of a given layer section (y direction). 
Then a ID Moiré pattern is produced in the orthogo- 
nal X direction, smoothly alternating between AA' , BA' 
and AB' stacking as shown in Fig. [l|a). The corre- 
sponding hopping amplitudes are related by Vaa' i^) — 
Vba'{x - L/3) = Vab'{x + L/3), where Vaa'{x) « 
{w/vf)[^ + 2 cos{2Trx/L)] using a single-harmonic ap- 
proximation [TH] (the interlayer coupling is w ~ t±/3 = 
0.11 eV, where t± is the interlayer hopping). 

To assist in interpreting the role of the different in- 
terlayer couplings, we deflne the functions Ax{x) — 
-{Vab'Ìx) + Vba'{x))/2 and Ay{x)^ = {Vab'{x) - 
Vba'{x))/2. Then Vab- = -A^ + Ay, Vba' = -A^-Ay, 



and it becomes clear that A^ , Ay act as off-diagonal vec- 
tor potentials. Taking Pauli matrices cr in the AB pseu- 
dospin space and r in the space of the two layers, we may 
recast Eq. ([T]) into 



H — vpcT ■ {—id — A) + vfVaa'Ti 



(2) 



where we have introduced the gauge potential A = 
{A^Ti, AyT2), which induces a precession of the layer in- 
dex as an electron moves in real space. This A is non- 
Abelian, since [A(r), A(r')] 7^ in general (see also Ref. 

This formulation highlights the different nature of 
the Vaa' coupling, that acts rather like a scalar potential 
(proportional to the unit matrix ao). 

This electron system has the characteristic property of 
developing fiat bands of spatially confined states at large 
L, whose formation is fully controlied by the effect of the 
gauge potential A. Computing the energy levels of the 
Hamiltonian ([2]), one observes that at large L the sys- 
tem develops two increasingly narrow subbands around 
zero energy of states confined between AB' and BA' re- 
gions, see Fig. |2] Their energy, for any given momentum 
kx and \ky\ < 3w/vp, oscillates towards zero, crossing 
it periodically as L increases (e.g. whenever wL/2TrvF 
is an integer if ky = 0, see inset on the right panel of 
Fig. [2]). Additionally, a second pair of fiat bands ap- 
pear at a finite energy, corresponding to states confined 
around AA' . Ali these bands become AA' -conUned and 
linearly dispersive in ky for \ky \ > 3w/vf, although they 
remain non-dispersive in the x direction. These features 
are strongly reminiscent of the Landau-level to snake- 
state transition in carbon nanotubes of large radius in a 
real perpendicular magnetic field "5^] , which have also an 
effectively modulated magnetic fìux. 

This confinement phenomenology may be understood 
from the effect of a confining potential created purely by 
the gauge field A. The equation for the eigenstates 
of H can be expressed after squaring the Hamiltonian 
(and disregarding for simplicity the scalar potential at 
this point) as 

{-d^ + id-A + 2iA-d + Al + Al~-a^Fxy)'Ì> = {e/vpf'f 

(3) 

where the field strength is conventionally defined in terms 
of the matrix-valued potential as = — 
d^A^ — i[Afj^,A^]. Given the invariance of H under 
the combined operation of charge conjugation and par- 
ity, the eigenstates can be chosen in the form ^'(r) — 
[0i(-r),(/)i(r),(/)5(-r),02(r)], for some (j)i^2- In the 
limit of zero transverse momentum ky, the combinations 
(f>±{'r) = 0i(r)±02(— r) decouple, and the above equation 
translates, at large L, into 

- vWiix) = -y±(x)<^±(x) -I- O (^) (4) 

wìthV^(x) = -{±e+Ax+Ay){±e+Ax-Ay) (SD]- Thisis 
the wave equation of a scalar mode with eigenvalue E = 
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FIG. 2. Left: Dispersion of the low-energy eigenstates of the 
Hamiltonian ([2| as a function of ky, for = Q and L ~ 
3700a, where a is the C-C distance. Note the zero-energy 
band (confined between AB' and BA') and its sateUite fiat 
band (confined around AA'). The inset covers a larger energy 
range. Right: Low-energy levels of the sheared bilayer as a 
function of the period L for kx — ky — . Note the two types 
of States, AA' -con&ned, which scale as e ~ 1/%/L (ali but the 
first, which scales as 1/L, see inset), and the AB'—BA' states, 
which cross zero energy when wL/vp = 27rn for integer n. 

under the infìuence of an e-dependent confìning potential 
V^{x), sketched in Fig. [l] e = eigenstates centered 
around AB' and BA' regions will arise whenever a level of 
such potential crosses E = 0. Such states will be peaked 
exactly at AB' and BA' , since the well has E — turning 
points at said regions. Moreover, a discrete set of i? = 
eigenstates centered around the AA' locai minimum will 
arise at energy e ^ 1/ ^/L. These two types of states are 
apparent in the numerical bandstructure plottcd on the 
right panel of Fig. [2] 

The above analysis in tcrms of Vcff relies crucially 
on the non-Abelian character of the gauge potential, 
[Ax,Ày] ^ 0. Without this property, the recurrence of 
zero-energy states as L increases would not appear. This 
may be appreciated from an alternative point of view. In 
order for a (normalizable) zero-energy state to exist, the 
operator Wg^o relating the wavefunction at x — and 
X = L, [0i(L),02(L)] = We=o[4>i{Q),4'2ÌO)], must have at 
least one eigenvalue of modulus one. Since at zero energy 
Eq. ^ leads to 

-^d^ ( ^2 ) " ^^''^ ^ " '^"^^^ {ti)' 
we have for ky = 

We=o = Pexp i^i J dx [Ax{x)ti ~ zAy(x)r2]| 

where "Pexp" denotes the path-ordered product of expo- 
nentials of differential line elements. [3T] One can check 
that this operator becomes unitary when wL/vp = 2TTn 
for integer n. This is the condition for the existence of 



normalizable zero-energy modes, in agreement with the 
numerical results. 

Low energy description of twisted bilayers. — At en- 
ergies below 1 eV, a twisted bilayer may be accurately 
modeled by Hamiltonian ([T]), where the shift AK in 
the relative position of the Dirac points in each layer 
Comes as a consequence of the rotation by the twist an- 
gle 9. If we take the originai position of the K points 
as K = (47r/3ao,0), the shift in each layer is given 
by ±AK/2 = (0,±A'sin(6l/2)). On the other hand, 
also fixes the size of the Moiré pattern unit celi, which 
grows as 6 decreases. More precisely, the Bravais super- 
lattice formed by the Moiré pattern has primitive vec- 
tors L± = L(v^/2,±l/2), where L = ao/2 sin(6'/2). 
This periodicity becomes exact on an atomic level when 
the rotation is commensurate and minimal, such that 
L = \/l + 3n + Sn^ao for some integer n > P^ . 

The interlayer coupling may be written in terms of 
a single periodic profile ^(r) = V{r + L+) = V{r + 
L_), in such a way that if we fix Vaa'{^) — V{r), then 
VAS'(r) - V{r + (L+ + L_)/3) and Vba'Ìv) = V{r - 
(L-|_ + L_)/3). A common procedure is to assume that 
the interlayer hopping is dominated by processes with 
momentum-transfer Qo = or equal to the reciprocai 
vectors Qi,2 = (±27r/%/3, 27r)/L [111 [20], so that ^(r) = 
(w/wp) exp(iQj • r). Coupling V is complex in this 
case, however, but we can stili carry out the procedure 
of the preceding section by defining Ai^ = —Re{VAB' + 
Vba')/2,A2, = lm{VAB' + VBA')/2,Aiy = lm{VAB' - 
Vba')/2 and A2y Re(VAS' - Vba')/2. We can then 
write the Hamiltonian for the twisted bilayer as 

H = VFcr-{-id-T3AK/2-A)+VF^ (5) 

with non-Abelian potentials A — {AixTi+A2xT2, AiyTi + 
A2yT2) and è = Re{VAA')Ti - Iui{Vaa')t2 

The mismatch AK of the Fermi points may be removed 
by carrying out a gauge transformation on the spinors, 
^I^ = exp((i/2)r3 AK • r)'ì, at the expense of introducing 
new potentials Vijir) = Vij(r)e~'^^ '". We finally get a 
modified expansion V^(r) — (w/vp) exp(iqj • r), with 

a star of three vectors q^. Note that |V^(r)| — |V^(r)| has 
conical singularities at the center of AB' /BA' regions. 

Two representative bandstructures obtained numeri- 
cally from the Hamiltonian ([s]) for different values of 6 
are plotted in Fig. |3] The first corresponds to an in- 
dex n = 20 and exhibits a lowest subband with vanish- 
ing energy at the two Dirac points originating from the 
graphene layers. As the angle 9 is decreased, the energy 
scale of the lowest subband is significantly lowered, until 
it becomes remar kably fiat for values of n around n = 31 
(0 « 1°), exhibiting zero Fermi velocity at the K point 
and a bandwidth that is more than 100 times smaller 
than the scale of the next subband. (Note however that 
this is not a topological zero mode in the sense of a stan- 
dard zero Landau level [3 3) , since the Atiyah-Singer index 
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FIG. 3. Left: Low-energy subbands of the Hamiltonian ([5| 
along the first Brillouin zone of the bilayer superlattice for 
n = 20 (dashed hnes) and n = 31 (full lines), for which a 
zero-energy band develops and the Fermi velocity at the K 
point vanishes. Right: Localization pattern (in logarithmic 
color scale, white is maximum) around A A' stacking of wave- 
functions on the zero-energy band for the first four values of 
n at which the Fermi velocity vanishes. 

[53] is zero). Lowering further, the lowest subband be- 
comes dispersive once more, before coUapsing again, and 
so on, showing a recurrent behavior as a function of the 
size L of the Moiré pattern [20] . 

For low values oi 6 {n > 31), the lowest-energy eigen- 
states show a strong confinement in the regions with AA' 
stacking, as shown in Fig. [Sj which is confirmed by 
atomistic tight binding calculations |l8j. This confine- 
ment is essentially controlied by the vector potential A, 
as the pattern of confinement remains unmodified when 
the scalar potential $ is ideally switched off in the model. 
The eigenstates obey now an equation similar to ([s]), but 
with Al+Al replaced by Al,^. + A^^ + A^y + A^y and Zee- 

man coupUng to F„j = dxAiyTi + dxA2yT2 - dyAi^Ti - 
dyA2xT2 + '2AixA2yTg, — 2A2xAiyT3. The contributions to 
the energy square of order ~ w'^ can be combined in the 
form {Ai-j-zL A2y)'^ + {A2xT Aiy)"^ . This function becomes 
zero only at the center of AA' stacking and at the center 
of either AB' or BA' stacking (depending on the eigen- 
values of az and T3). This degeneracy is broken by the 
derivative terms in F^y, which tend to confine at points 
where the gradients of Vab' and Vba' become higher. 
These functions become flatter at the regions of AB' and 
BA' stacking, respectively, and are more steep at the cen- 
ter of AA' stacking, explaining the effect exerted by the 
vector potential to confine in the latter region. 

We note that the first instance at which the low- 
est subband becomes fiat has a simple interpretation as 
the situation where the analogue of the magnetic length 
Ib ~ y/vpLjw starts to fit in the bilayer supercell of 
size L. One can actually check that, ai n ~ 31, the re- 
sult of computing the flux integrai ip = J (PrF^y leads 
to values (p ~ <ì>o'''2, <ì'o(cos(7r/6) ri — sin(7r/6) T2), and 



— <I>o(cos(7r/6) Ti + sin(7r/6) T2) for supercells rotated by 
2tt/3 in the twisted bilayer, with $0 = Stt (in units 
h ^ 1). This corresponds to the flux quantum rotated 
in the SU(2) flavor space. Unlike for that first instance, 
higher values of n giving rise to a fiat lowest-energy sub- 
band do depend on the strength of the Vaa' coupling 
[35] . However, the essential spatial confinement prop- 
erties of the corresponding lowest-energy eigenstates do 
not. They remain confined around AA' stacking. They 
also acquire higher angular momentum components and 
become increasingly ring-shaped for higher values of n 
(see Fig. [3]), as expected for the excited states of a 2D 
potential well centered around AA' stacking. 

Experimental measures of the low-energy electronic 
properties of twisted bilayers have been reported in par- 
ticular in Ref. It has been found that, at a certain 
value 6^1°, the renormalized Fermi velocity near the 
K point of the twisted bilayer becomes so small that the 
picture based on Dirac quasiparticles breaks down. This 
Comes together with the observation of a clear pattern of 
spatial confinement in the locai density of states, which 
adopts the form of a triangular charge density wave fol- 
lowing the modulation of the Moiré pattern. These fea- 
tures are fully consistent with the confinement of the low- 
energy eigenstates in the regions of AA' stacking due to 
the action of the gauge potential, which provides a strong 
confinement mechanism according to the preceding dis- 
cussion. This single-particle mechanism will cooperate 
with the additional many-body effects that may also con- 
tribute to the modulation of the charge in the system. 

Conclusion. — We have shown that the Moiré-like 
modulation of the interlayer hopping in graphene bilay- 
ers leads to a very rich phenomenology, which can be de- 
scribed in terms of effective non-Abelian gauge potentials 
in the low-energy electronic theory. We have shown that 
any additional terms arising from the stacking modula- 
tion, such as non-Abelian scalar potentials, do not qual- 
itatively modify the low energy electronic structure. In 
the case of sheared bilayers with quasi-lD Moiré patterns, 
the gauge potential is equivalent to a confining potential 
that leads to low-energy charge accumulation along ID 
strips. The effect of the non-Abelian gauge potential in 
rotationally faulted bilayers is also to develop a character- 
istic spatial pattern of confinement, and the formation of 
dispersionless bands for discrete value of the Moiré peri- 
ods. We conclude these two effects are the characteristic 
signature of Moiré-induced non-Abelian gauge potentials 
in graphene bilayers. 

The emergence of these types of gauge fields is generic 
to Systems of coupled Dirac equations, and the analysis 
presented bere can be extended to multilayered systems 
with SU(N) gauge groups. One may also furthcrmore 
envision the possibility of tuning the non-Abelian fields 
caused by stacking by applying generic strain fields to 
Moiré bilayers. These will not only give rise to Abelian 
fields as in monolayers, but also to small modifications 
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of the stacking non-Abelian fields [53] , whose interplay is 
known to produce a rich phenomenology (15) . 
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